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Driving magnetic skyrmions with microwave fields 
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We show theoretically by numerically solving the Landau-Lifshitz-Gilbert equation with a classical 
spin model on a two-dimensional system that both magnetic skyrmions and skyrmion lattices can be 
moved with microwave magnetic fields. The mechanism is enabled by breaking the axial symmetry 
of the skyrmion, for example through application of a static in-plane external field. The net velocity 
of the skyrmion depends on the frequency and amplitude of the microwave fields as well as the 
strength of the in-plane field. The maximum velocity is found where the frequency of the microwave 
coincides with the resonance frequency of the breathing mode of the skyrmions. 
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Skyrmions, topologically stable magnetization textures 
with particle-like properties, have recently attracted 
great attention due to their potential use in future 
spintronic devices [5]. The manipulation of skyrmions 
is of great importance and interest: skyrmions can 
be driven using spin-polarized current EH3, magnetic 
or electric field gradients mm , temperature gradi¬ 
ents [12UI4] and magnons USHUj. Microwaves, on the 
other hand, have been broadly used in studying vari¬ 
ous magnetic phenomena, such as the ferromagnetic reso¬ 
nance (FMR) and spin wave excitations in skyrmion crys¬ 
tals [T8H2T] . However, the possibility of creating transla¬ 
tional motion of skyrmions has not been explored in these 
experiments [ISM]. In this Rapid Communication, we 
show that both a single skyrmion and a skyrmion lattice 
can be moved by microwave fields if the axial symme¬ 
try of skyrmions is slightly broken by a static in-plane 
external field. 

We employ skyrmions stabilized by the Dzyaloshinskii- 
Moriya Interaction (DMI) [22 113] ■ More precisely, the 
bulk DMI is considered so that a chiral skyrmion (vortex¬ 
like) rather than a hedgehog (radial) skyrmion configura¬ 
tion emerges [2 US US] • We start with a classical Heisen¬ 
berg model on a two-dimensional regular square lattice 
with nearest-neighbor symmetric exchange interaction, 
the bulk-type DMI, and the Zeeman field PH CEB [25] . in 
addition, a time-dependent magnetic field h(i) is applied 
in the +z-direction. Accordingly, the system’s Hamilto¬ 
nian can be written as 

H ~ J E(ij) 1,1 • m, + •> ' hi, x nij] 

- Ei H( H + h(t)) • Hi.;, (1) 

where (i,j) represents a unique pair of lattice sites i 
and j, m, is the unit vector of the magnetic moment 
Hi = —fr/Si with Si being the atomic spin and y(> 0) 
the gyromagnetic ratio, and J is the symmetric exchange 
energy constant. In the case of bulk DMI, the DMI vec¬ 
tor D ij can be written as D , 7 = Dr^, where D is the 
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DMI constant and r ^ is the unit vector between and 
Sj. We use the DMI value with D/J = 0.18, which re¬ 
sults in the spiral period A ~ 2nJa/D ~ 25 nm for a 
typical lattice constant a = 0.5 nm [S]. We consider two 
nonzero components for the static external field H: an 
in-plane component H y and a perpendicular component 
H Zl i.e., H = (0 ,H y ,H z ). A nonzero H z is essential for 
stabilizing the skyrmion crystal [15]. 
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FIG. 1. (a) Skyrmion configuration in the presence of an in¬ 
plane field H y = 0.006 with D = 0.18, J = 1 and H z = 0.02. 
(b) The symmetric topological charge density distribution of 
a skyrmion when H y = 0. (c) The corresponding topological 
charge density for the skyrmion shown in (a) when H y > 0. 
(d) The skyrmion lattice with 12 skyrmions in a sample of 
size N = 174 x 150 sites, with an in-plane field H y = 0.004. 


The spin dynamics at lattice site i is governed by the 
Landau-Lifshitz-Gilbert (LLG) equation, 

dm t dm, 

, )f = x H 0 ff + am, x (2) 

where a is the Gilbert damping and H c ff is the effec¬ 
tive field that is computed as H cff = — (l/|/i i |)d'H/dm i . 
The Hamiltonian (|T|) associated with the LLG equa¬ 
tion ([2]) can be understood as a finite-difference-based 
micromagnetic model. Therefore, our simulation results 
are reproducible by setting the saturation magnetization 
M s = h"/S/a 3 , exchange constant A = J/2a and DMI 
constant for continuum form D a = —D/a 2 (correspond¬ 
ing to the energy density £dmi = D a m- (V x m)) in micro- 
magnetic simulation packages such as OOMMF |27| . We 
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TABLE I. Unit conversion table for J = 1 meV, 5=1 and 
a = 0.5 nm. 


Distance x 

x = a 

— 0.5 nm 

Time t 

t = hs/j 

« 0.66 ps 

Velocity v 

v = Ja/(hS) 

« 7.59 x 10 2 m/s 

Frequency ui 

Cj = J/{hS) 

« 1.52 x 10 3 GHz 

Magnetic field H 

H = J/(hyS) 

« 8.63 T 


have carried out simulations with and without dipolar 
interactions, and the results are qualitatively the same. 
We report results without dipolar interactions for clarity 
of the model assumptions. 

A two-dimensional system of size N = 160 x 160 
sites with periodic boundary conditions is selected to 
study the dynamics of a single skyrmion, Fig. 0 a ) , and 
N = 174 x 150 sites for the skyrmion lattice, as shown 
in Fig. [l|d). We have chosen J = H = r y = S = a= l 
as simulation parameters [18], 25J, therefore, the coeffi¬ 
cients to convert the external field H , time t, frequency 
u> and velocity v to SI units are H = J/fvyS, t = hS/J , 
Cj = J/hS and v = Ja/hS , respectively. Table I shows 
the expressions, and particular values for the case of 
J = 1 meV, 5=1 and a = 0.5 nm. We use simulation 
units throughout the paper. The perpendicular compo¬ 
nent H z is fixed as H z = 0.02 which corresponds to 0.173 
T for S = 1 and J = 1 meV. We use Gilbert damp¬ 
ing a = 0.02 for all simulations except for the magnetic 
spectra shown in Fig. [2] where a = 0.04 is chosen. For 
the single skyrmion dynamics, we apply the absorbing 
boundary conditions for damping [29j by setting a = 1.0 
for the 20 spins at the edges of a simulated domain. 

The configuration of a skyrmion in the presence of an 
in-plane field H y = 0.006 is shown in Fig.[lj(a). It is found 
that the radial symmetry is broken. Indeed, as shown in 
Fig [ljc) , the corresponding distribution of the topolog¬ 
ical charge density q(x,y) = (l/47r)m • ( d x m x d y m) is 
asymmetric. However, the total topological charge of a 
single skyrmion remains constant Q = J qdxdy = — 1. As 
a comparison, Fig. 0 b ) shows the topological charge den¬ 
sity q for a skyrmion with radial symmetry when H y = 0. 
Similar to the vortex [301 , the distortion of the skyrmion 
is along the cc-axis when an external field is applied in 
the y-direction. 

The excitation of internal modes depends on the static 
external field H z as well as the frequency and direc¬ 
tion of microwaves mi eh- The typical excited modes 
are the clockwise/counterclockwise rotation and breath¬ 
ing modes j3BH20]. To study how the in-plane ap¬ 
plied field H y affects the excitation mode of a skyrmion, 
we calculate the magnetic absorption spectrum of the 
skyrmion. After applying a sinc-function field pulse h = 
hosinc(uJot) = ho sm(u>ot) / (uiot) to the stable skyrmion 
state we record the spatially averaged magnetization evo¬ 
lution and from that we compute the dynamic suscepti¬ 
bility x via a Fourier transformation [321 55] . For in¬ 
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FIG. 2. Imaginary parts of the Hy-dependent dynamical sus¬ 
ceptibility as a function of frequency for (a) a single 
skyrmion, and (b) a skyrmion lattice. The spectra are ob¬ 
tained by applying a sine-field pulse h = hosinc(uoot) to the 
system with ho = 1 x 10 -5 and uoo = 0.17T in the 2 -axis, the 
magnetization dynamics is recorded every dt = 5 for 8000 
steps. 


stance, the component Xzz is computed using m z when 
the pulse is parallel to the 2 -axis. 

Figure [2] shows the imaginary part of the dynami¬ 
cal susceptibility Xzz for a single skyrmion (a) and a 
skyrmion lattice (b); each calculated for different in-plane 
fields H y . We see in (b) for the skyrmion lattice and 
H y = 0 that the mode with frequency ui ~ 0.0246 is dom¬ 
inant. This is the so-called breathing mode [18l EE]. The 
resonance angular frequency ui ~ 0.0246 in simulation 
units corresponds to the frequency / = oj/2tt « 5.95 GHz 
(using Cj from Table [I]) . The breathing mode frequency 
decreases slightly as the in-plane field H y increases. A 
second peak emerges with increasing H y , and the fre¬ 
quency of the new mode is w ~ 0.0278 when H y = 0.006. 
Similar to the skyrmion lattice, a breathing mode with 
ui ss 0.0168 is found for a single skyrmion, as shown in 
Fig. [2](a). In general, the resonance frequency of a sin¬ 
gle skyrmion is lower than that for the skyrmion lattice. 
As for the skyrmion lattice, a new mode with frequency 
around 0.02 emerges for the single skyrmion case as H y 
is increased. This second mode is the uniform mode with 
frequency to = jH where H ~ (H y + H z ) 1 ^ 2 is the am¬ 
plitude of the external field. 

In the presence of an in-plane applied field H y , the 
skyrmion is deformed, as shown in Figs. [l](a) and[l](c). 
Therefore, instead of the geometric center we measure the 
so-called guiding center |51] R = {X,Y) of a skyrmion: 
X = J xqdxdy/ f qdxdy and Y = J yqdxdy/ f qdxdy , 
where q is the topological charge density. For a symmet¬ 
ric skyrmion, the guiding center is the same as its geo¬ 
metric center. In the rest of this work, we consider the 
scenario that a linearly polarized microwave is applied in 
the z-direction, i.e., h(t) = ho sin(wot)e z , where ho and 
uio are the amplitude and frequency of the microwave, 
respectively. 

Figure [3](a) shows the displacement of the guiding 
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FIG. 3. (a) The displacements of the guiding center (X, Y) 

for a single skyrmion. The simulation parameters are loo = 
0.017, ho = 2 x 10 -4 and the in-plane field H y = 0.004. (b-d) 
The velocities v x and v y as functions of (b) the microwave 
frequency wo, (c) the microwave amplitude ho, and (d) the 
in-plane field H y . The fixed simulation parameters are the 
same as in (a). 


center for a single skyrmion with u>q = 0.017 and in¬ 
plane field H y = 0.004. The microwave amplitude is 
h 0 = 2 x 10 , which corresponds to 1.73 mT for J = 1 

meV and 5 = 1. It can be seen that the x-component of 
the guiding center, A', changes significantly as a function 
of time, while the displacement of Y is relatively small. 
Figure [3^b) plots the frequency-dependent skyrmion ve¬ 
locity and shows that a single skyrmion has maximum 
velocity when uiq = 0.017, which is the breathing mode 
resonance frequency. Therefore, exciting the breathing 
mode can move the skyrmion effectively in the presence 
of an in-plane field H y . Using the conversions presented 
in Table I, the maximum velocity is v x ~ 2.8 cm/s. While 
v x is positive for frequency Wo = 0.017, it is negative for 
wo = 0.023, where the former corresponds to the breath¬ 
ing mode, and the latter to coherent rotation. The re¬ 
lation between the velocities and the amplitude of the 
microwave is shown in Fig. [3jc) , and the dependence of 
v x on ho is nonlinear: v x oc which is proportional to 
the power of the microwaves. 

Figure §d) describes the relation between the 
skyrmion velocity and the in-plane field H y . The velocity 
is zero if H y = 0, which is expected due to the symme¬ 
try of the skyrmion. The velocity of the skyrmion also 
depends on the direction of the in-plane field H y : the ve¬ 
locity is reversed when the direction of the in-plane field 
is reversed. Similarly, a change in the sign of the DMI 
constant will also reverse the sign of the velocity, which 
is different from the case of driving skyrmions with spin- 
polarized currents, where the sign of perpendicular veloc¬ 
ity (with respect to the current direction) of the skyrmion 
motion is related to the sign of topological charge rather 
than the DMI constant sign. 

To understand why the skyrmion moves in the pres¬ 
ence of an in-plane field, we split the magnetization 



FIG. 4. The total spatial force density fi = m s ■ [d;m s x (m x 
H c ff)] for (a) H y = 0, and (b) H y = 4 x 10 —3 , where we have 
used m s = (m). The microwave frequency is u>o = 0.017. 


unit vector m into a slow part m s and a fast part n, 
i.e., m = m s + n, where the slow part represents the 
equilibrium profile of the skyrmion while the fast part 
is responsible for the excited spin wave mode [121135]. 
In the continuum approximation, the effective field is 
H e g = AV 2 m — Z)V x m + H + h(f), where A = 2 A/M s 
and D = 2 D a /M s . In the presence of the microwaves 
with specific frequency, bound spin-wave modes are ex¬ 
cited, and thus we expect (n) = 0 and (n) = 0 due to 
the microwave synchronization, where T = 2ir /wq is the 
period of microwaves and the notation (/) = ( f)(t) = 
T _1 • J f(t')dt' represents the time average of func¬ 
tion f(t) over a single period T. Furthermore, one ob¬ 
tains (m s x h) « 0 and (n x rn s ) ss 0 since m s is the 
slow part. Therefore, by averaging the LLG equation ©> 
over a period T we arrive at 

(m s ) = -7( m x H off ) + a(m s x rh s ), (3) 

where (n x ri) = 0 is used since basically n is a sine 
or cosine function in time. We then consider the possible 
translational motion of the skyrmion such that m s (r, t) = 
m s (r — v s t), i.e., rh s = — (v s • V)m s , where the skyrmion 
velocity v s = dR/dt is assumed to be a constant. If the 
skyrmion moves slowly, i.e., v s T <C L (L is the typical 
skyrmion size), we have (m s ) « — (v 3 • V)m s (see Ap¬ 
pendix |A| where m s = (m s ). Similarly, (m s x m s ) ss 
m s x (m s ), and thus Eq. ([3]) can be rewritten as 

(v s • V)m s = 7 (m x H off ) + am s x (v s • V)m s . (4) 

Following Thiele’s approach in describing the motion of 
magnetic textures [36) . we replace the dots in J m s • 
((3;m s x • • • )dxdy by Eq. Q to obtain mmm 

G x v s + Vv s = F, (5) 

where i = x,y and G = 47rQe z . The tensor = arjij is 
the damping tensor in which rjij = J(difh s -djih s )dxdy = 
dijT) is the shape factor of the skyrmion and rj is close to 
47T [35]. The force F is given by 

Fi = -7 J m s • [<9im s x (m x H off )] dxdy. ( 6 ) 
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Figure |4^a) and (b) depict the total spatial force den¬ 
sity for H y = 0 and H y = 4 x 10 —3 , respectively, where 
we have used m s = (m). The force density is sym¬ 
metric if H y = 0 and thus the total force F is zero. 
However, when H y is nonzero the force distribution is 
asymmetric which results in the skyrmion motion due to 
the nonzero net force. For small damping a <C 1, we 
have v x « F y /(AnQ). The total force calculated with 
parameters ujq = 0.017, H y = 0.004 and ho = 2 x 10 ~ 4 
is F y = —4.7 x 10 —4 , therefore, the established veloc¬ 
ity is v x = 3.7 x 10 -5 , which fits the simulation result 
(~ 3.8 x 10~ 5 ) well. Similarly, for to = 0.023 using Eq. M 
we obtain Fy = 5.7 x 10 -5 and find v x ~ —4.5 x 10^"° 
from Eq. (15]); in agreement with the simulation results 
(the minimum of v x is —4.7 x 1CU 6 ). 

It is of interest to circumstantiate the contributions of 
the total force F. In Appendix |B] we show that there are 
three nontrivial terms 

(m x H off ) « (n x [AV 2 n-OTxn + h (*)]). (7) 

The exchange term nx AV 2 n corresponds to magnon cur¬ 
rents nans]. Compared to the skyrmion motion induced 
by the temperature gradient, where the magnon cur¬ 
rent is generated by the temperature gradient, here the 
magnon current originates from the external microwave 
fields. Another difference is that in our case the contribu¬ 
tions from DMI and microwave fields are also significant. 
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FIG. 5. (a) The velocities v x and v y as a function of an 

in-plane field H y for the skyrmion lattice at frequency to o = 
0.023. (b) The velocities as a function of frequency too with 
H y = 4 x 10~ 3 . The microwave amplitude is ho = 2 x 10 -4 . 


We repeat the velocity study for the skyrmion lattice. 
Fig. ]5](a) plots the velocities v x and v y of the skyrmion 
lattice as functions of the in-plane external field H y . The 
dependencies are similar to the single skyrmion case. The 
frequency-dependent velocities v x and v y are shown in 
Fig. §b). As for the single skyrmion case, the velocity 
peak coincides with the dominant dynamical susceptibil¬ 
ity peak in Fig. |2jb) . 

In closing, we briefly comment on the importance of 
symmetry breaking in driving the skyrmions. The driv¬ 
ing force originates from the microwave field, which is 
periodic in time and averages to zero. The symmetry¬ 
breaking field converts the periodic microwave field into 
a net force and thus moves skyrmions effectively. A re¬ 
lated field with periodic driving forces is that of ratchet¬ 


like transport phenomena 1371 - 139] . where the net motion 
is obtained by breaking the spatial symmetry |37j or tem¬ 
poral symmetry [ 35] . We also note that preliminary sim¬ 
ulation results suggest that for magnon-driven skyrmions 
manEi the introduction of a symmetry-breaking in-plane 
field affects the skyrmion’s motion and changes the Hall 
angle significantly. 

In summary, we have studied the skyrmion dynamics 
driven by microwaves in the presence of an in-plane ex¬ 
ternal field. We found that both a single skyrmion and 
a skyrmion lattice can be moved by a linearly polarized 
microwave held if the axial symmetry of skyrmions is 
slightly broken. These results suggest a novel method 
for skyrmion manipulation using microwaves fields. 
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of Southampton, in the completion of this work. 


Appendix A 

Assume that a well-behaved function f(x,t) = f(x — 
vt) describes the dynamics of a soliton where v is a con¬ 
stant. As we can see, / satisfies / = —vf. For given 
time T, if vT <C L where L is the typical size of the soli¬ 
ton (for example, L could be the domain wall width for 
a magnetic domain wall), we can find that 

(/>('0) = fdt= ^[f{x - vT) - f(x)} 

*-vf'(x-vT/2), (Al) 

where we have used the Taylor series for fix) and f(x — 
vT): 

fix) « fix - vT/2) + f\x - vT/2)vT/2 (A2) 

f{x - vT) « fix - vT/2) - fix - vT/2)vT/2.iA3) 

Similarly, we can see that f = (/)(0) ~ fix — vT/ 2) and 
thus we have (/)(0) ~ —vfo- This relation actually holds 
for arbitrary t 

if) « -«/■ (A4) 


Appendix B 

By using the effective field explicitly and noticing that 
m = m s + n, the term (m x H c ff) can be splited into 
four parts 

(m x H off ) = <Ti + T 2 + T 3 + T 4 >, (Bl) 

where Tf = n x [4V 2 n — l)V x n + h(f)] is shown in 
Eq Q, T 2 = m s x (AV 2 m s - Z)V x m s + H), T 3 = 
n x (AV 2 m s — DV x m s + H) and X 4 = m s x [dV 2 n — 
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ZJV x n+h(f )]. We expect = 0 since m s represents the 
equilibrium state of the skyrmion. For the slow skyrmion 
motion, replacing m s by m s and noticing that (n) = 0 , 


we obtain (T 3 ) ~ (n) x (AV 2 m s — DV x m s ) = 0 and 
(T 4 ) ss m s x (AV 2 n — ZJV x n + h(£)} = 0. In this slow 
motion approxmation, the fast part n can be computed 
as n « m (m). 
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